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Diagnosing the burst influence upon accretion in the clocked burster GS 1826-238 
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ABSTRACT 

Type-I X-ray bursts on the surface of a neutron star are a unique probe to the accretion 
in X-ray binary systems. However, we know little about the feedback of the burst emission 
upon accretion. Hard X-ray shortages and enhancements of the persistent emission at soft 
X-rays have been observed. To put these findings in context with the aim of understanding 
the possible mechanism underneath, we investigated 68 bursts seen by RXTE from the clocked 
burster GS 1826-238. We diagnosed jointly the burst influence at both soft and hard X-rays, and 
found that the observations can be described as the CompTT model with variable normalization, 
electron temperature and optical depth. Putting these results in a scenario of coronal Compton 
cooling via the burst emission would lead to a shortage of the cooling power, which may suggest 
that additional consideration like the influence of the burst on the corona formation should be 
accounted for as well. 

Subject headings: stars: coronae - stars: neutron - X-rays: binaries-X-rays: bursts 


1. Introduction 


A low-mass X-ray binary (LMXB) consists of a compact star (i.e., a black hole or a neutron star) and 
a donor. The donor usually fills its Roche lobe, and thus transfers mass to th e compact star via releasing the 
gravitational potential energy as persistent X-rays (e.g., iLewin et al.ll2006h . Its spectrum can be described 
as a combination of a thermal and a non-thermal (power-law like) components. The thermal component is 
considered as the emission from an optically thick accretion disk, a nd the powerdaw component is believed 


to originate from hot plasmas above or in side the accretion disk (IWhite 


speculated as being the "corona" (see e.g. iFrank et al.l l2002l: Done et al 


1988 


2OO7I: 


Mitsuda 


19891). which is 


Zhandl2013h . In case of 


having a neutron star (NS) as the compact object, the mass accumulated on the NS surface may end up with 
a violent thermonuclear explosion, in which the integrated flux co uld reach up to typ i cally 10^^-10^" erg s. 


These intense X-ray flares are called type-I X-ray bursts (see e.g. Lewin et al.lll993l:lGallowav et al.ll2008h . 
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Being a short intense shower of soft X-rays, type-I X-ray bursts can be used as a probe of the LMXBs 
(low mass X-ray binaries) accretion via investigating the evolution of the persistent emission along with 
the runaway of the thermonuclear explosion. The bursts are expected to have a blackbody spectrum wit h 
slight deviations that depend on the burst luminosity (for details, see, e.g., ISuleimanov et al.ll2011al. 120121) . 
To extract the burst spectrum, one usually assumes that the p ersistent emission is u nchanged during th e 
burst and one subtracts it off as part of the background (see e.g. lKuulkers et al.ll2003t) . IWorpel et al.l (l2013n . 
however, proposed that the persistent emission intensity is variable; thus involving a multiplicative factor /a 
to the persistent emission for fitting the spectrum during the burst (assuming an unchanged spectral shape, 
which we shall refer hereafter as the varying persistent flux model). They found that in photospheric radius 
expansion bursts (PRE bursts) the best-fit value of the /a is significantly great er than 1, while it seem s that 
the /a values are ~ 1 when the persistent emission is near the Eddington limit, lin’t Zand et al.l (l2013h found 
spectral excesses from a blackbody spectrum at both low and high energies in SAX J1808.4-3658, which 
can be eliminated by using the varying persistent flux model. The value of /a seems to be proportional 
to the int ensity of bursts , while this correlation may reverse at higher burst luminosity in the hard state of 
EMXBs (IJi et al.ll2014bh. In addition, enhanced persistent emission was also reported during a superburst 
from 4U 1636-536 (IKeek et al.ll2014ah . This shows the significant influence of the burst emission upon the 
accretion at soft X-rays (~ 2-10 keV), although the underlying mechanism is still not understood. Possible 
explanations may be Povnting-Rober tson effect (IWalkenll992h or some other reprocessing of burst photons 
by the accretion disk and the corona din’t Zand et al.ll2013h . In the hard X-rays (e.g., above 30keV), on the 
contrary, the pers istent emission during the burst is found to have a clear shortage, indicating cooling of a 
corona (see, e.g.. iMaccarone et al.ll2003al: IChen et al.ll2012l: |ji et al.ll2014ah and hence dramatic changes of 
the persistent spe ctrum, since the bu rst emission dominates below 30 keV. This issue was partially addressed 
in previous work ( Keek et ^l2014a l. They performed a detailed time-resolved spectral analysis for this rare 
superburst event and found evidences in the decay phase of an enhancement of the overall persistent flux and 
a reduction of the corona temperature, which could be consistent with a hard X-ray shortage if the persistent 
spectrum continues up to tens of keV. Since the superburst event is relatively rare, we explore this issue 
in more details in this paper by analysing the atoll source NS EMXB GS 1826-238, which shows more 
frequent X-ray bursts. 

GS 1826-238 is a persistent atoll source, wh i ch does not exhibit si gnificant spectral evolution of its 
persistent emission until recently (lNakahirall2014l: iGallowav et al.ll2008h . Assuming a dis tance of 6 kpc. 


this s ource always stays in the so-called low hai'd state with a luminosity of 5%-9% Ledd (IGallowav et al. 


20081). Its X-ray emission in a wide energy range of 3-200 keV can be well desc ribed by a pure the rmal 
Comptonization model with a characteristic hot electron temperature, ~ 20keV dCocchi et ahlboiol) . GS 
1826-238 is the fa mous "clocked burster" because of the recu rrence of the bursts being relatively stable over 
the span of years (lUbertini et al.ll 19991: iGalloway et al.ll2004l) . Also the burst shapes remain quite constant, 
and the burst peak fluxes have a rather narrow distribution (see, e.g.. Table 1). The mean and standard 
deviation of the burst peak fluxes are 24.6 and 2.8, respectively. These properties are quite helpful for 
reducing the possible systematic errors when one stacks bursts together. The most distinguishing feature for 
GS 1826-238 is that it is one of a few atoll XRBs with a hai'd X-ray shortage being significantly detected 
during bursts (IJi et al.ll2014al) . Since for GS 1826-238 the thermal emission from either accretion disk or 
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NS surface contributes little at energies above 30 keV, such a hard X-ray shortage provides the most clean 
probe to the change of the non-thermal emission, which likely helps to restrict the spectral model at soft 
X-rays with respect to the pre-burst emission. 

In this paper we derive the evolution of the persistent emission along the burst at soft X-rays, and then 
diagnose the burst influence on accretion in a broad energy band, using the knowledge of the previously 
detected hard X-ray shortage. 


2. Observations and data analysis 


We analyzed all Rossi X-ray Timing Explorer (RXTE) Proportional Counter Array (PCA) pointing 
observations and searched for burstsli] The bursts exhibiting data losses were excluded. We also ruled out 
several bursts for which the Standard 2 mode data were not available . 43 bursts wer e previously chosen for 
studying the hard X-ray shortage in the type I bursts of this source (IJi et al.ll2014ah . Since in this analysis 
the selection criteria are a bit looser, i.e., bursts with incomplete persistent emission (150 s before and 250 s 
after the burst peak) are also considered, 25 more bursts come into the overall sample. All of the 68 bursts 
selected for our analysis are shown in Table 1. 


In this paper we analyzed the persistent emission using Standard 2 mode data , while the Event mode 
data were used for studying the properties of bursts because of their high time resolution. The software 
environment was HEAsoft ver 6.15, including the commands "saextrct" and "seextract" to extract light 
curves and spectra and "pcabackest" to estimate the background. The spectral response files were generated 
by the latest pcarsp VI 1.7.1, in which the drifting PCU gain values over the RXTE mission lifetime have 
been taken into account. We have verified that the results are not sensitive to selection of the RXTE epoch. 
The spectral analysis was performed by using Xspec 12.8 and its Python wrapper PyXspec. During spectral 
fits, the used energy band was 3-20 keV, above which the effective area of PCU drops off rapidly. To avoid 
the possible bias in estimation of parameters because of the non-Gaussian distribution in some channels with 
small numbers of counts, a Churazov weighj^ was used in the spectral analysis. 


We per formed time-res olved analysis on each burst, by adopting the same spectral fitting procedure as 
described in iji et all (l2014bl) . In brief, we performed time-resolved analysis for the individual bursts around 
their peaks. To improve the statistics, we selected variable time intervals so that the total counts in each time 
bin exceeds 5000 photons. As a result, the bin size has to gradually increase toward the burst tail because 
of the lower count rate. The varying persistent flux model, i.e., wabs * {bbodyrad + fa x persistent flux) 
was employed to fit the spectra in the 3-20 keV band, in which the bbodyrad represented the i ntrinsic burst 


comp onent, while /a x persistent flux represented the variable persistent emission (see, e.g.. lWorpel et al. 


201-31 . for details). The "wabs" component was used to describe the interstellar absorption; the hydrogen 


'To search for the bursts, we produced 1-s light curves, then select the observations in which the maximum count rate is 500 
cts/s larger than the average, and finally checked them manually. 

^see https://asd.gsfc.nasa.gov/XSPECwiki/low_count_spectra 
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column density was fixed to 0.3 x 10^^ atoms cm“^ dCocchi et al.ll2010h . Our results are independent of 
the photoelectric absorption model due to the fact that its major influence happen at energies below 3 keV. 
The varying persistent flux model was chosen because the short burst exposure does not allow measuring all 
parameters of the relatively weak persistent emission (see Table 1). Here the persistent flux was estimated 
by fitting the spectrum before each burst with a phenomenological model, such as wabs * (powerlaw+gauss) 
or wabs *{bbodyrad+powerlaw+gauss), where the Gaussian component was always set at 6.4 keV with an 
arbitrary width of 0.5 keV. More physical models, e.g., CompTT, would lead to large parameter errors due 
to RXTE/PCA’s limited bandpass, but was addressed in the subsequent simulations. An average reduced-x^ 
of ~ 1.3 (45 dof) shows that it is appropriate to describe phenomenologically the persistent emission with 
the adopted model. 


3. Results 


3.1. /a evolution and correlation with burst flux 


As mentioned above, we performed time-resolved analysis for individual bursts to obtain the curves of 
the burst flux and /a factor. For individual bursts, however, the possible change of the /a parameter is not 
very significant due to relatively large error bars (see Figure 1 in Lli et al.l (l2014bh '). In order to improve the 
statistics, we investigated the average evolution of /a by summing up all the bursts. In practice, we took the 
peak flux of each burst as a time reference (t=0 s) and averaged over bins of 3 seconds. Figure [T] shows that 
at the rising phase of bursts, the /a value rapidly increases (O to A). Then, along with a further increasing of 
the average burst flux, the /a value decreases and reaches a value of around 1 at time of the burst flux peak 
(noted as point B in Figure[Tll. Finally, during the averaged burst decay, /a gradually returns to a value (point 
C) comparable to that at point A, and then smoothly decreases. We note that if we divide the bursts into 
two groups along time and then calculate /a curves individually to check their consistency, the x^ test leads 
to a re duced-51.83(46 doH, which indica tes that the slight changes of the shape of bursts over several 
years (iGallowav et al.l l2004l: lGallowayll2012h have little influence on our results. For comparison. Fig. [U 
also presents the aver aged hard X-ra y GO-50 keV) lightcurve during bursts observed by RXTE/PCA using 
the Event mode data (I.Ti et al.ll2014ah . □ The hard X-ray shortage is significant around the burst peak and 
gradually recovers to its initial level. 


To illustrate the /a evolution against flux, we present these data also in the Figure |2] The correlation is 
significantly different between the cases when the burst flux is larger and smaller than ~ 10“^ erg s“^ cm“^. 
The turnover is smooth instead of a sharp peak. The Pearson correlation coefficient is 0.87 with a p-value 
of 2.4 X 10 for the former, but -0.95 with a p-value of 2.6 x 10 ^ for the latter. This result strongly 


confirms what was reported previously: /a is highly dependent on the burst luminosity (IWorpel et al.ll2013 


^We note that compared with the hard X-ray lightcurve, 25 more bursts are used to produce the soft X-ray and the /a lightcurves 
with smaller error bars (i.e., upper panel in Figure 1). The soft X-ray and the /a lightcurves are changed little if we use identical 
selection of bursts. 
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in’t Zand et alJl2013l : iKeek et alJl2014al). Although a quite similar /a trend was reported previously in the 
hard state of 4U 1608-522 (IJi et alJl2014hh. the fact that no significant shortage of the persistent hard X-ray 
emission is observed in this source (IJi et alJl2014ch . prevents us from a joint diagnosis of the /a evolution 
over a broad energy band. 


As shown in 


■Ti et al.l (l2014al') for GS 1826-238, the hard X-ray shortages for the persistent emission 


are still significant around points A and C of the bursts, suggesting a cooled corona. However, at that 
time, the /a curve at soft X-rays show large enhancements of the persistent emission (/a ~ 1.8). Around 
the burst peak, both the hard X-ray flux and the /a present a similar trend of anti-correlation with the 
burst flux, which likely suggesfs a common cause. An e nhancemenf in fa can be fhe re sulf of, e.g., fhe 
infrinsic spec frum fhaf deviates from a blackbody slighflv (ISuleimanov ef al.ll2011aL l2012h or fhe reflection 
componenfs (IBallanfyne ef al.ll2004l. 120051: iKeek et al.ll2014bh . which might lead to a biased estimate of /a 
when the "bbodyrad" model is used. 

In what follows we diagnose these possible mechanisms in the light of the possible burst influence on 
the persistent emission, using the results of /a as der ived at yft X-r ays in this paper and of the persistent 
emission shortage reported previously at hard X-rays (IJi et al.ll2014all . 


3.2. Diagnosing the burst influence 


Here we adopt the "burstatmo" model to describe the intrinsic spectra of bursts instead of the simple 
bbodyrad fsee. ISuleimanov et al.ll2011aL 1201211 . The "burstatmo" model represents the atmospheres of neu¬ 
tron stars and emergent spectra in a plane-parallel, hydrostatic, and local thermodynamic equilibrium (LTE) 
approximation. Hence, we performed again the above-described analysis, but now using a slightly different 
model: wabs * {burstatmo+fa x persistent flux). There are three parameters in the burstatmo model, i.e., 
the radius and luminosity of the neutron star, and a normalization. We can not simultaneously constrain 
these parameters due to the large error bars and chose to fix the less impacting parameter-the radius of 
the neutron star-at 10 km. If the increasing value of /a up to ~ 1.8 originates from the deviation of the 
intrinsic burst spectrum from a blackbody, after using "burstatmo" model, fa should be around 1. However, 
we found that the resulting fa is very similar to the case when bbodyrad model was used (i.e., similar to 
the Figure [T]but with slightly bigger error bars). In addition, the "burstatmo" model does not improve the 
goodness of fit. The reduced is 1-10 (18 dof), which is comparable with 1.08/18(dof)using the model 
wabs*{bbodyrad+fa x persistent flux). In order to further rule out the possibility that the /a > 1 originates 
from the bias of fitting due to the spectral subtle deviations between the intrinsic spectra (i.e., burstatmo) 
and phenomenological bbodyrad model, we used the burstatmo model to produce faked spectra with as¬ 
sumed burst luminosities 0.5 and 1, respectively, and then fitted them by varying persistent flux model using 
bbodyrad model. The resulting best-fit values of fa are 1.31 and -4.01, respectively, which are inconsistent 
with the observed values (i.e., 1.6 and 1.1 with error bars of ~ 0.1) shown in Figured] 

Alternatively, the reflection du ring bursts produced by the a c cretion disk repres ents viable possibilities 
for influencing the persistent flux (IBallantvne et al.ll2004l. l2005l: iKeek et al.ll2014bh . These were detected 
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durin g superbursts din’t Zand et alJl2010l : iKeek et alJl2014ah or hinted at in normal bursts (iGallowav et al. 
20081). For individual bursts, the statistics are not good enough to search for possible reflection components. 
Therefore, we stacked the residuals derived from the spectral analysis using the varying persistent flux 
model. If the reflection component does exist, one would expect some structure in the residuals. The 
analyses were performed at three instances (t=0 s, 20 s, 60 s) with a time bin of ~ 1 s for the different phases 
of the (average) burst. We found, however, no additional structures suggestive of emission lines or edges. 
The excess around 6.4 keV is hardly visible with a significance level below 2 cr. 

The h ard X-ray shortage during bursts in GS 1826-238 (and other LMXBs) is indicative of cooling of 
the corona (IMaccarone et al.ll2003al: lji et al.ll2014al) . The result reported from a superburst in 4U 1636-536 
also showed evidence in the burst tail for changes of the electron temperature in the corona when the 
blackbody normalization was fixed dKeek et al.ll2014al) . We investigated whether the /a variability shown in 
Figure m corresponds to any corona cooling by performing simulations. We set free the electron temperature 
Tg and produced some faked spectra, and fitted them using the varying persistent flux model to look into 
a possible influence on /a factor. We used the model CompTT to describe the persi stent flux; the mode l 
parameters (i.e., ro=1.02, re=21.5keV and r=3.81) were adopted from the report of ICocchi et al.l (l201(]l) . 
We note that a slightly different model, e.g., cutoffpl or comptb, barely affects our analysis. Assuming the 
spectrum during the bursts is wabs * (bbodyrad+CompTT) with variable Tg, we produced faked spectra by 
using the Xspec command/ake/t, where the bbodyrad components are set at typical values T=2keV and 
norm= 100. Then we used the varying persistent flux model to fit these faked spectra, i.e., assuming the 
spectrum remains constant in shape but varies in normalization. The resulting value /a, as illustrated in 
Figure [3l clearly indicates that a lower temperature Tg tends to result in a smaller /a. 


Considering that the 3 additional independent variables in the coronal model CompTT, i.e., the normal¬ 
ization (A), the optical depth (r), and the temperature of the seed photons (To), are also expected to vary 
under the shower of soft X-rays from bursts, we further investigated if the /a evolution can be connected to 
variations in these parameters. For simplicity, we first kept r and Tq unchanged during bursts, but set free the 
variables N and Tg; see below for a discussion of variations in the other parameters. We performed simula¬ 
tions by producing faked spectra based on the model wabs * {bbodyrad+N x CompTT{Tg)), in which there 
are two free variables: N and Tg. We again used the varying persistent flux model to fit these faked spectra, 
like done before. The normalization of CompTT is fixed to 4.78 x 10“^, for which the inferred counts rate is 
consistent with the hard X-ray (30-50 keV) observations, i.e., 1.69 cts s“^ pcu“'. The assumed temperature 
and normalization in the bbodyrad model are 2keV and 100, respectively, which corresponds to the values 
around the burst peaks; note that the simulations are insensitive to the assumed blackbody parameters. The 
result is illustrated in the left panel of Figure IH The colorbar in that figure shows the value of /a, whereas 
the black lines show the contours where /a equals to 0, 1 and 2, respectively. We note that in Figure |4] the 
/a is just a result of using varying persistent flux model, and does not reveal physical meanings, i.e., the 
enhanced (suppressed) accretion rate when the /a is larger (smaller) than 1. The result suggests that larger N 
and Tg lead to a larger /a value. This is not surprising since both of the N and /a are a normalization; and Tg 
will affect the efficiency of the up-Compton scattering. Thus, the scenario of corona cooling via bursts can 
be sketched as follows. At the beginning of the burst, the normalization (x-axis, N) increases rapidly, which 
























-7- 


phenomenologically accounts for the /a > 1 (0->A, point O represents the initial values before the bursts, 
i.e., Te = 2\.5keV and N=\). This process may be accompanied by a moderate decrease of the Tg. Then, Tg 
drops sharply and dominates the burst influence by largely reducing /a (A->B and B->C). In the cooling 
tails the process might be a rough inverse of the 0->A process, but with a dominance of corona heating over 
than the cooling0 We note that the burst flux and the /a value at point A are comparable with these at point 
C, while its hard X-rays seems to be slightly stronger, which may imply a little larger Tg at point A. 


A s mentioned ab ove, the hard X-ray shortage in the bursts of GS 1826-238 is observed to be rather 
robust ( Ji et alJboiTa ). This provides a unique opportunity for testing the cooling scenario speculated above, 
via checking if the observed hard X-ray shortage is consistent with the persistent spectral model derived from 
the varying persistent flux model. For a given persistent model and assumed parameters, one can deduce 
the hard X-ray count rates by convolving it with PCA’s response matrix, and then compare to the results 
using observations directly. Therefore, the observed hai'd X-ray shortage can impose a strong constraint in 
the parameter space of the persistent spectral model. The hard X-ray count rate at 30-50 keV before the 
bursts is ~ 1.69 ± 0.1 cts s“^ pcu“', and is reduced to a minimum ~ 0.8 ± 0.1 cts s“' pcu“^ around the burst 
peak. Assuming the coronal spectrum is A x CompTT(Te), these hard X-ray shortages correspond to the 
upper and lower cross hatches in the left panel of Figure |4] respectively, between which is the constrained 
parameter space of N and Tg. We found that for points A and C the parameter spaces having hard X-rays 
between ~ 0.8 ±0.1 cts s“^ pcu“* and ~ 1.69 ± 0.1 cts s“* pcu“^ (i.e., between the two hatches) do not 
overlap the region in which the /a is ~ 1.8. This implies that the CompTT model with variable N and Tg can 
not reconcile the observations at both hard and soft X-rays. 


To solve such an incompatibility, one may consider additionally the other two variables in CompTT, 
i.e., the temperature of seed photons To and the optical depth r. We found that the temperature of seed 
photons influences little both the hard X-ray shortage and the /a evolution, probably because Tq (which is ~ 
1 keV) is comparable to the burst temperature ( ~ 1-2 keV). r , however, has a great eff ect on the persisten t 
spectrum, especially at the hard X-rays. When r decreases from an initial value of 3.81 dCocchi et al.ll2010n 
to ~ 2, the parameter region inferred from the hard X-rays shortage (at 30-50 keV) is compatible with that 
estimated from /a evolution (at 3-20 keV). This result is shown in the right panel of Figure IH and hints at a 
variable optical depth r as a consequence of the burst. 


4. Discussion and Summary 


Feedback of type -I X-ray bursts onto the accretion material surrounding the compact star in LMXBs has 
been predic t ed earlier ( Walker et al.lll989 : WalkeJl992 : Ballantvne et ahlboOd . 2005 : Inogamov & Sunvaev 


1 999l. l20 1 d : iKluzniakI|20 1 3l) . However, due to the generally short exposure of bursts, it is hai'd to test them 


in individual events. Using RXTE observations of GS 1826-238 spanning seventeen years, we were able to 


"'Note that Figure|4]is just a schematic diagram and the positions of the points (A/C, B) are not precise; they only representing 
the values of /a ~ 1.8, and ~ 1, respectively. 







































study in detail the overall c hanges of the persisten t emission while bursting. This study can be done due to its 
constant burst shape (e.g.. [Thompson et al.ll2005n and dominance of the coronal emission in th e hard state, 
and, therefore, sign ificant non-thermal radiation above ~ 30keV (e.g.. iGallowav et al.ll2004ll . Following 
Worpel et all (1201311 . we employed the varying persistent flux model to fit the spectra in the energy band 3- 


20 keV during bursts. The /a is correlated with the burst flux up to ~ 10 * erg s * cm^ (which corresponds to 
~ 0.8 X 10^^ erg s“^ assuming a distance of 6 kpc and a gravitational redshift of 1.3), while ab ove this value 


it is anti-correlated. This confirms the previous report using the hard state of 4U 160 8-522 (Woroel et al. 


20131: lji et al.ll2014bh . In 4U 1608-522, however, no hard X-ray shortages were found (IJi et al. 


2014ch 




We found that, although the varying persistent flux model is successful in fitting the persistent spectrum 
during bursts, especially in the 3-20 keV energy band, more physical mechanisms are required to account 
for the contemporary hard X-ray shortage. We considered a variety of possibilities among which only 
the inverse Compton cooling of the corona seems reasonable. We found that an increase of the coronal 
model compTT normalization (N) by a factor of 2-3 can account for the /a enhancement with respect to 
the pre-burst emission. This increased N might resu lt from the in creasing seed photons provided by type-I 
X-ray bursts, or the Poynting-Robertson drag effect (IWalkenll992h . We note that the luminosity of the burst 
overwhelms that of the pre-burst accretion disk by a factor of more than 10, while the CompTT normalization 
only increases by a factor of 2-3. This suggests that the interaction between the accretion dis k and the corona 
is more efficient, which would be the case if the c orona "sandwic hes" the disk ( Whitelll988 ). but not when it 
is spheri cal in the inner side of the truncated disk ( Mi tsudall 1989 1. This result is in agreement with the recent 
report of izhan 3 (12^ . who found that the spectra of dipping (i.e., high-inclination) LMXBs, regardless of 
their spectral states, exhibit stronger comptonized emission, suggesting that coronae may have disk-like or 
oblate shapes. 

Around the point B in Figure [T] /a has a trend similar to that of the strength of the hard X-ray shortage. 
This hints at a sufficiently cooled corona when the burst emission is intense. Figure |4]shows that the electron 
temperature Tg may account for the /a valley observed around the burst peak. Such a T- reduction/recovery 
feature was observed as well in the decay phase of the superburst from 4U 1636-536 (IKeek et al.ll2014ah . 
which also points to corona cooling under a shower of soft X-rays. We find, however, that to reconcile the 
hard X-ray shortage with the /a evolution at soft X-rays, one has to take into account a smaller optical depth 
r. The Comptonised spectrum can be used to estimate the cooling power, i.e., the power of losing energy 
of the corona under the shower of burst photons. Here we estimated the cooling power as the difference 
between the incident flux and the emergent flux after the Compton process. The emergent flux is estimated 
by using Xspec command "flux" with the parameters specified belo\\@, and the incident flux is inferred by 

^^Max{T C) 

analytic estimation, i.e., the emerg ent energy of photons is ~ ’ times larger than the incident 

energy (IRvbicki & Lightmanlll979h . If we assume that around the burst peaks (point B) Tg ~ 13keV and 


^We note that the lack of a hard X-ray shortage might he caused by e.g., the cont amination of bursts having high color tempera¬ 
ture, or the lack of significant non-thermal emission. For details, see iji et all ( l2014cl) . 

®The "flux" command can be used to calculate the flux of the specified model, parameters and energy ranges. For details, see 
https;//heasarc.gsfc.nasa.gov/xanadu/xspec/manual/XSflux.html 
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r ~ 2, the inferred cooling power is 6.8 x ergs in the energy band 1-100 keV, which is only ~ 

8% of the cooling power before bursts. The uncertainties of the cooling power is difficult to be constrained. 
If we artificially assume fhe error bars of and r are ~ 3 keV and 1 af 90% confidence levels, respectively, 
and fhey are independenf, fhe 5% lower and upper boundaries of fhe inferred disfribufion of fhe cooling 
power are 1.5% and 24.0% of fhe cooling power before bursfs, respectively. This c ooling power furns ouf fo 
be nof sufficienf fo accounf for fhe decreased femperafure of elecfrons in fhe corona ( Keek ef al .boUah . Thus 
we speculafe fhaf a possible solufion mighf be fo fake info accounf an additional mechanism, i.e., fhe corona 
formation, which may be influenced as well by fhe bursf. This may be reminiscenf of fhe XRB oufbursfs 
in fheir specfral fransifion from hard sfafe fo soff sfafe, where fhe corona is evenfually diminished along 
wifh fhe increasing accretion rafe. If one assumes fhaf fhe changes of fhe elecfron femperafure are fhe resulf 
of fhe balance befween fhe Compfon cooling and heating mechanism, a suppressed heating process mighf 
be a solufion fo fhe above confradicfion. However, fhe heafing mechanism currenfly is still poorly known. 
In evaporation models, bofh fhe disk and corona are individually powered by fhe release of gravifafional 
energy associafed wifh fhe accr efion of mafter affecfed fhrough viscous sfresses, which depends sfrongly 
on fhe viscosify pai'amefer (a) (IMeven 1 994l: iLiul l20 1 3h . Al fernafively , in magnetic reconnecfion models, 
fhe heafing power is ~ where Va is fhe Alfven speed (ILiul 120021) . Therefore, if af any poinf in time 
fhe sysfem is in approximafe equilibrium, i.e., fhe heafing power equals fo fhe cooling power, fhe viscosify 
paramefer a or fhe Alfven speed Va is expecfed fo change dramatically during bursfs. We caufion fhaf 
fhis scenario mighf be oversimplified. In oufbursf fransifions, which is also affribufed as fhe cooling and 
reheafing of fhe corona, fhe low sfafe fo high sfafe specfral fransifion af fhe inifial rise of oufbursfs occurs 
af a luminosify 5 times g reafer fhan fhaf of fhe fransifion from high sfafe fo low sfafe af fhe declining sfage 
(IMaccarone ef al.ll2003bl) . which is known as "hysferesis effecf" . This means fhaf fhe seed phofons’ flux 
mighf be nof fhe only factor fhaf can influence fh e cooling and re heafing of fhe corona, and some ofher 
unknown processes should be faken info accounf dRegelmanl 120141) . We speculafe fhaf during fransifions 
from fhe hard sfafe fo fhe soff sfafe in oufbursfs of LMXBs, fhe corona is cooled by fhe increasing soff seed 
phof ons emiffed from fhe accre fion disk, which in general leads fo an increased value of r (for defails, see, 
e.g.. lNaravanlll995l : lEsinlll997h . This seems inconsisfenf wifh whaf we infer for GS 1826-238. As shown in 


fhe leff panel of Figure |4j fhe main inconsisfency fhaf resulfs from /a and fhe hard X-ray shorfage is af fhe 
poinfs A and C. A smaller optical depfh of fhe corona has fo be infroduced fo accounf for fhe difference of fhe 
persisfenf emission modifications befween fhe soff and fhe hard X-rays. We speculafe fhaf fhis ’discrepancy’ 
mighf arise from fhe differenf opfical depfhs of fhe corona for fhe seed phofons having differenf origins. 
Considering fhaf if fhe corona is geomefrically fhick above fhe accrefion disk, buf exisfs only wifhin a small 
radial radius, fhe opfical depfh for fhe seed phofons from fhe accrefion disk would be larger for fhese phofons 
from fype-I X b ursfs. In addi tion, fhe seed phofons from fype-I bursfs would lead fo significanf Poynfing- 
Roberfson drag dWalkenl 19921) . This m ay imply fhaf during bursfs fhe inner part of fhe corona vanishes and 
becomes fhinner in fhe radial direction (IWalkenll992h . leading fo a smaller opfical depfh. 
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Table 1: The columns denote the observation ID, the time when type-I bursts occurred, persistent flux and 
unabsorbed peak flux of bursts at 3-20 keV. The persistent flux, peak flux of bursts and their errors are 


obtained by fitting models with Xspec command "cflux". 


Observation ID 

Modified Julian Day 

Persistent Flux 
(10“’ergs-' cm“’) 

Peak Flux 
(10“^ erg s“^ cm“^) 

Observation ID 

Modified Julian Day 

Persistent Flux 
(10“’ erg s“' cm“2) 

Peak Flux 
(10“^erg cm“^) 

30054-04-02-000 

50971.70 

0.996 ±0.010 

28.8 ± 0.9 

90043-01-01-01 

50971.70 

1.507 ± 0.005 

25.5 ±0.8 

30054-04-02-01 

50971.23 

1.012 ±0.004 

27.7 ± 0.9 

90043-01-01-02 

50971.23 

1.581 ± 0.007 

26.5 ± 0.8 

30054-04-02-02 

50972.18 

1.126 ±0.007 

28.1 ± 0.9 

90043-01-01-020 

50972.18 

1.577 ± 0.013 

26.0 ± 0.8 

30054-04-03-02 

50976.42 

1.014 ±0.004 

27.5 ± 0.7 

90043-01-02-00 

50976.42 

1.617 ± 0.009 

24.7 ± 0.9 

30060-03-01-01 

50988.83 

0.991 ±0.006 

27.1 ± 0.7 

90043-01-02-001 

50988.83 

1.578 ± 0.026 

25.6 ±0.9 

50035-01-01-01 

51724.88 

1.309 ±0.007 

24.6 ± 0.7 

90043-01-02-001 

51724.88 

1.578 ± 0.026 

24.2 ± 0.8 

50035-01-02-00 

51725.71 

1.176 ±0.019 

27.3 ± 0.7 

90043-01-02-01 

51725.71 

1.636 ± 0.007 

24.7 ±0.8 

50035-01-02-00 

51725.88 

1.176 ±0.019 

26.8 ± 0.7 

90043-01-02-01 

51725.88 

1.636 ± 0.007 

25.1 ±0.8 

50035-01-02-02 

51726.72 

1.364 ±0.009 

26.0 ± 0.7 

90043-01-03-000 

51726.72 

1.620 ±0.021 

26.5 ±0.8 

50035-01-02-04 

51728.77 

1.359 ±0.014 

26.0 ± 0.7 

90043-01-03-01 

51728.77 

1.551 ± O.OII 

25.7 ±0.8 

50035-01-03-00 

51811.75 

1.324 ±0.016 

26.5 ± 0.8 

90043-01-03-01 

51811.75 

1.551 ± O.OII 

24.7 ±0.8 

50035-01-03-08 

51813.49 

1.361 ±0.012 

24.9 ± 0.7 

91017-01-01-000 

51813.49 

1.537 ± 0.028 

23.5 ±0.8 

50035-01-03-09 

51814.01 

1.417 ±0.017 

26.4 ± 0.8 

91017-01-01-05 

51814.01 

1.600 ±0.015 

25.7 ±0.8 

50035-01-03-10 

51814.35 

1.397 ±0.007 

26.0 ± 0.8 

91017-01-02-01 

51814.35 

1.580 ±0.011 

23.8 ±0.8 

50035-01-03-120 

51813.15 

1.368 ±0.011 

26.1 ± 0.7 

91017-01-02-010 

51813.15 

1.543 ± 0.006 

25.9 ±0.8 

70044-01-01-00 

52484.57 

1.578 ±0.017 

27.6 ± 0.8 

91017-01-02-010 

52484.57 

1.543 ± 0.006 

24.8 ±0.8 

70044-01-01-000 

52484.42 

1.717 ±0.038 

24.7 ± 0.8 

91017-01-02-03 

52484.42 

1.587 ± 0.009 

26.6 ±0.8 

70044-01-01-02 

52485.01 

1.615 ±0.024 

26.4 ± 0.7 

91017-01-02-09 

52485.01 

1.583 ± 0.015 

25.7 ±0.8 

70044-01-02-00 

53205.20 

1.580 ±0.013 

24.3 ± 0.8 

91017-01-02-11 

53205.20 

1.557 ± 0.045 

25.2 ±0.8 

70044-01-02-01 

53206.06 

1.514 ±0.012 

26.3 ± 0.8 

92031-01-01-000 

53206.06 

1.257 ± 0.008 

20.2 ± 0.6 

70044-01-04-00 

53966.37 

1.654 ±0.010 

24.5 ± 0.8 

92031-01-01-000 

53966.37 

1.257 ± 0.008 

20.8 ± 0.7 

80048-01-01-00 

52736.53 

1.417 ±0.009 

23.6 ± 0.7 

92031-01-01-001 

52736.53 

1.227 ± 0.014 

17.8 ±0.6 

80048-01-01-010 

52736.67 

1.410 ±0.006 

24.4 ± 0.8 

92031-01-01-010 

52736.67 

1.251 ± 0.005 

17.9 ±0.5 

80048-01-01-010 

52736.80 

1.410 ±0.006 

24.0 ± 0.6 

92031-01-01-010 

52736.80 

1.251 ± 0.005 

18.7 ±0.5 

80048-01-01-04 

52738.48 

1.348 ±0.007 

24.5 ± 0.7 

92031-01-01-011 

52738.48 

1.283 ± 0.024 

20.5 ± 0.7 

80048-01-01-07 

52738.61 

1.397 ±0.014 

23.6 ± 0.7 

92031-01-01-011 

52738.61 

1.283 ± 0.024 

19.3 ±0.7 

80048-01-01-07 

52738.75 

1.397 ±0.014 

32.4 ± 0.9 

92031-01-01-011 

52738.75 

1.283 ± 0.024 

20.2 ± 0.7 

80049-01-01-00 

52820.56 

1.507 ±0.009 

24.1 ± 0.5 

92031-01-02-00 

52820.56 

1.282 ± 0.013 

20.2 ± 0.6 

80049-01-02-00 

52944.40 

1.425 ±0.010 

26.4 ± 0.8 

92031-01-02-000 

52944.40 

1.231 ± 0.007 

19.8 ±0.7 

80049-01-03-02 

52834.39 

1.515 ±0.010 

25.1 ± 0.8 

92031-01-02-000 

52834.39 

1.231 ± 0.007 

19.6 ±0.6 

80049-01-03-03 

52835.31 

1.523 ±0.013 

26.9 ± 0.7 

92703-01-04-02 

52835.31 

1.561 ± 0.009 

24.4 ± 0.8 
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Fig. 1.— The evolution of the averaged flux, value of /a (upper panel) and hard X-rays at 30-50 keV (bottom 
panel). The red horizontal line represents that the count rate at 30-50 keV equals to 1.69 cts pcu~^. The 
green vertical lines represent that the time equals to -8, 0, 20 s, respectively. 
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Fig. 2.— The burst flux vs. the value of /a. The burst flux is estimated as A'KaT'^N, where a is the Stefan- 
Boltzmann constant, T and N are the best-fit parameters of blackbody model "bbodyrad", respectively. The 
dash line represents the flux equals to 0.8 x 10^^ erg s“^ 
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Fig. 3.— The expected /a vs. the electron temperature T^. The dashed lines show the initial parameters 
before bursts, i.e., re=21.5 and /a=l. The red lines represent the characteristic error bars when stacking the 
burst spectra with exposure time 0.5 s, 1 s and 2 s, respectively. 



Fig. 4.— The resulting /a by assuming the r = 3.81 and 2 for the left and right panel, respectively. The 
colorbar shows the intensity of the /a and the black lines shows the contours where the /a equals to 0, 1 and 
2, respectively. The two cross hatches e nclose the poss ible region of N and the Tg, which are deduced from 
the hard X-rays shortages in Figure 2 in Lli et al.l (l2014ah . The green lines show the initial state before bursts, 
/a =1 and re=21.5 keV. 
































































































































































